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Abstract 

We have measured the branching fraction £>, longitudinal polarization fraction f L , and CP 
asymmetry coefficients A and S for B° — > p + p~ decays with the Belle detector at the KEKB 
e+e" collider using 2531b" 1 of data. We obtain B = 22.8 ± 3.8 (stat) ±|;| (syst) 1 x 10~ 6 , f L = 

0.941 t lt (stat) ± 0.030 (syst), A = 0.00 ± 0.30 (stat) ± 0.09 (syst), and S = 0.08 ± 0.41 (stat) ± 
0.09 (syst). These values are used to constrain the Cabibbo-Kobayashi-Maskawa phase 2 ; the 
solution consistent with the Standard Model is 4> 2 = (88 ± 17)° or 59°<0 2 <115° at 90% C.L. 

PACS numbers: 13.25.Hw, 12.15.Hh, 11.30.Er 
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One of the main goals of the e + e~ "5-factories" is to determine whether the Cabibbo- 
Kobayashi-Maskawa |l| mixing matrix with three quark generations is unitary; failure to 
satisfy this criterion would indicate new physics. Unitarity imposes six independent con- 
straints upon the matrix elements, one of which is V* b V ud + V* b V cd + V^ b V td = 0. Plotting 
this relationship in the complex plane yields a triangle, and unitarity is tested by measuring 
the internal angles (denoted 1; 2 , 3 ) to check whether they sum to 180°. The angle 2 is 
the phase difference between VXV td and —V* b V ud and is measured via b— >u decays such as 
B ^n + n~, p ± 7r =F , and p + p~ [2|. Of these, B°^>p + p~ gives the most precise value as the 
contribution from a possible loop amplitude (with a different weak phase) is smallest. The 
size of the loop amplitude is constrained by the upper limit on B{B°^p°p ) 0. 

One determines 2 by measuring the At distributions of B°B° events, where At is the dif- 
ference between the decay time of the signal B° (B°) and that of the opposite-side B° (B°). 
For B°/B° — > p + p" decays, these distributions have interference terms of opposite sign 
proportional to e~' A *' // ' r s[^4cos(Am At) + S sin(Am At) ], where Am is the B°-B° mass 
difference and A, S are functions of 2 . Here we present a measurement of the B°-^p + p~ 
branching fraction B, longitudinal polarization fraction f L , and coefficients A and S, using 
253 fb _1 of data recorded by the Belle experiment p| at KEKB 0j. 

Candidate B° — > p + p~ , p^ — >7r ± 7r° decays are selected by requiring two oppositely charged 
tracks satisfying p T > 0.10 GeV/c, dr < 0.2 cm, and \dz\ < 4.0 cm, where p T is the 
momentum transverse to the beam axis, and dr and dz are the radial and longitudinal 
distances, respectively, between the track and the beam crossing point. The tracks are fitted 
to a common vertex. We require that tracks be identified as pions based on information 
from a time-of-flight system, an aerogel Cerenkov counter system, and the central tracker p| . 
The resulting identification efficiency is about 89%, and the kaon misidentification rate is 
about 10%. Tracks are rejected if they satisfy an electron identification criterion based on 
information from an electromagnetic calorimeter (ECL). 

The ^ candidates are combined with ir° candidates reconstructed from 7 pairs hav- 
ing M 77 in the range 117.8-150.2 MeV/c 2 (±3cx in m^o resolution). We require E > 
50 (90) MeV in the ECL barrel (endcap), which subtends 32°-129° (17°-32° and 129°-150°) 
with respect to the beam axis. To identify p ± — > ^ir decays, we require that M^ ±it0 be 
in the range 0.62-0.92 GeV/c 2 (±2T in the M w±n0 distribution). To reduce combinatorial 
background, the 7r°'s must have p > 0.35 GeV/c in the e + e~ center-of-mass (CM) frame, 
and p candidates must satisfy —0.80 < cos# ± < 0.98, where 9± is the angle between the 
direction of the 7r° from the p and the negative of the B° momentum in the p^ rest frame. 

To identify B° — > p + p~ decays, we calculate variables M bc = J E£ eam — p 2 B and AE = E B — 
E heam , where E B and p B are the reconstructed energy and momentum of the B candidate, 
and -Ebeam is the beam energy, all evaluated in the CM frame. The AE distribution has a 
tail on the lower side due to incomplete containment of the electromagnetic shower in the 
ECL. We define a signal region M^G (5.27, 5.29) GeV/c 2 and AE e (-0.12, 0.08) GeV. 

We determine whether a B° or B° evolved and decayed to p + p~ by tagging the b flavor of 
the non-signal (opposite-side) B decay in the event. This is done using a tagging algorithm [f| 
that categorizes charged leptons, kaons, and A's found in the event. The algorithm returns 
two parameters: q, which equals +1 ( — 1) when the opposite-side B is most-likely a B° {B°); 
and r, which indicates the tag quality as determined from Monte Carlo (MC) simulation and 
varies from r = for no flavor discrimination to r = 1 for unambiguous flavor assignment. 

The dominant background is e + e~ qq (q = u, d, s, c) production. We discriminate 
against this using event topology: e + e~— > qq events tend to be jet-like in the CM frame, 
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while e + e~ — > BB tends to be spherical. To quantify sphericity, we calculate 16 modified 
Fox- Wolfram moments and combine them into a Fisher discriminant 0]. We calculate a 
probability density function (PDF) for this discriminant and multiply it by a PDF for 
cos#g, where 9 B is the polar angle in the CM frame between the B direction and the 
beam axis. BB events have a 1 — cos 2 9 B distribution while qq events tend to be uniform in 
cos#g. The PDFs for signal and qq are obtained from MC simulation and a sideband [M bc e 
(5.21, 5.26) GeV/c 2 ], respectively. These PDFs are used to calculate a signal likelihood C s 
and qq likelihood C qq , and we require that 71 = CJ (C s + C q g) be above a threshold. As 
the tagging parameter r also discriminates against qq events, we divide the data into six r 
intervals (denoted £=1 — 6) and determine the 1Z threshold separately for each. 

The overall efficiency (from MC simulation) is (3.19 ± 0.02)%. This value corresponds 
to f L = l; the change in efficiency (+5.0%) for f L equal to its central value measured below 
is taken as a systematic error. The fraction of events having multiple candidates is 9.5%; 
most of these arise from fake 7r°'s combining with good tracks, and thus we choose the best 
candidate based on \M —m n0 \. In MC simulation this correctly identifies the B° — > p + p~ 
decay about 90% of the time. A small fraction of signal decays (5.7% for longitudinal 
polarization) have > 1 7r ± daughters incorrectly identified but pass all selection criteria; 
these are referred to as "self-cross-feed" (SCF) events. Their vertex positions (and hence At 
values) are smeared. 

We determine the signal yield using two unbinned maximum likelihood (ML) fits. We 
first fit the M hc -AE distribution in the wide range M bc e (5.21, 5.29) GeV/c 2 and AE e 
(—0.20, 0.30) GeV to obtain the B° — > (p + p~ + nonresonant) yield N^ pp+nr y, we then fit 
the ¥^±^0 distribution of events in the M hc -AE signal region to obtain the nonresonant 
p ± 7r T 7T° + 7r ± 7T =F 7r 7r° fraction. 

For the first fit we include PDFs for signal p + p~ and 6— >c, b—+u, and qq backgrounds. 
The PDFs for signal and b-^u are two-dimensional distributions obtained from MC simu- 
lation; the PDF for b— >c is the product of a threshold ("ARGUS" (H) function for M bc and 
a quadratic polynomial for AE, also obtained from MC simulation. The PDF for qq is an 
ARGUS function for M bc and a linear function for AE; the latter's slope depends on the tag 
quality bin i. All qq shape parameters are floated in the fit. The b— >u background is dom- 
inated by B— > (p7r, a x 7r, a l p) decays; as their contributions are small, their normalization 
is fixed to that from MC simulation. For B + — > afir and B — > a x p modes, the branching 
fractions (unmeasured) used in the simulations are 3 x 10~ 5 and 2 x 10~ 5 , respectively; we 
vary these by ±50% and ±100%, respectively, to obtain the systematic error due to these 
estimates. The result of the fit is N^ pp+nr ^ = 207 t 2 g events. Figure [T] shows the final event 
sample and projections of the fit. 

For the subsequent fit, we require that events be in the M hc -AE signal region and fit 
w o in the wide range 0.30—1.80 GeV/c 2 . One p candidate is required to satisfy £ 
(0.62, 0.92) GeV/c 2 ; the mass of the other p candidate is then fit. We include additional 
PDFs for nonresonant B-^pirn and B^-kh-kh decays; these are taken from MC simulation 
assuming three- and four-body phase space distributions. However, the fit result for unim 
is <Cl% and thus we set this fraction to zero. The PDFs for p + p~ and b — > u are also 
taken from MC simulation. The PDFs for b — > c and qq are combined and taken from the 
sideband M bc G (5.22, 5.26) GeV/c 2 ; we check with MC simulation that the shapes of these 
backgrounds and their ratio in the sideband region are close to those in the signal region. 
We impose the constraint that the fraction of (p + p~± pirir) events in the M^±^ range 0.62- 
0.92 GeV/ c 2 equals that obtained from the M hc -AE fit; there is then only one free parameter. 
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FIG. 1: (a) M bc for events with AE £ (-0.10, 0.06) GeV. (b) AE for M bc G (5.27, 5.29) GeV/c 2 . 
The curves show fit projections: the dashed curve is p + p~ + p7T7r, the dotted curve is qq, the 
dot-dashed curve is b — >c, the small solid curve is b^u, and the large solid curve is the total. 




M m (GeV/c<) 
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FIG. 2: (a) M ± for events in the M hc -AE signal region that satisfy M (not fit) € 
(0.62, 0.92) GeV/c 2 . (b) Sum of cos#_|_ distributions for events in the signal region that satisfy 
M ± (both) £ (0.62, 0.92) GeV/c 2 . The curves show fit projections: the dashed curve is p + p" , 
the dot-dashed curve is p7rir, the dotted curve is qq + (£>— >c) + (6— >n), and the solid curve is the 
total. 



The fit obtains f fJ7T7T = f pn7T /U' PP + f P nn) = (6-3 ±6.7)%, and thus N pp = (1 - f p7T7T )N [pp+nr) = 
194 ±32, where the error is statistical and obtained from a "toy" MC study (since the errors 
on f p7T7T and N^ pp+n ^ are correlated). This value agrees well with the p + p~ yield obtained 
from the M^o fit (141 events) multiplied by the ratio of acceptances (1.33). Figure |2fa) 
shows the data and projections of the fit. 

The branching fraction is N J(e ■ e n ■ N -^), where N is the number of B° — > p + p~ 
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candidates, N is the number oiBB pairs produced [(274.8 ±3.1) x 10 6 ], e is the acceptance 
and event selection efficiency obtained from MC simulation, and e n is a correction factor 
for the ir^ identification requirement to account for small differences between data and the 
simulation (0.969±0.012). The result is B = (22.8±3.8) x 10" 6 , where the error is statistical. 

There are eleven main sources of systematic error. These are typically evaluated by 
varying the relevant parameter(s) by la and noting the change in B. The sources are: track 
reconstruction efficiency (1.2% per track); ir° efficiency (4% per 7r°); calibration factors 
(obtained from a large B + — > D°p + — > A' + 7r~7r°p + sample) used to correct the signal M bc - 
AE PDF to better match the data; the M hc -AE shapes for b— >c; the fraction and M hc -AE 
shapes for 6— >«; the AE range fit; statistics of the MC simulation used to calculate e; the 
dependence of s upon the polarization; uncertainties in e n and N B ^; and the qq suppression 
requirement. Combining these in quadrature gives a total systematic error of +10.1% and 
-11.6%. Thus, 



22.8 ± 3.8 (stat) (syst) x 10~ e . (1) 



To determine the longitudinal polarization fraction f L , we perform an unbinned 
ML fit to the 9 + , 9_ helicity angle distribution. This distribution is proportional to 

4f L cos 2 9, cos 2 9_ + (1 — f L ) sin 2 6 + sin 2 9_ . In the fit, this PDF is multiplied by an accep- 
tance function determined from MC simulation. The acceptance is modeled as the product 
A(cos8 + ) ■ A(cos8_), where A is a polynomial. 

We fit events in the M hc -AE signal region that satisfy M T±lfl e (0.62, 0.92) GeV/c 2 . 
We include PDFs for signal, pirn, and b — ► c, b — > u, and qq backgrounds. The PDFs for 
b^c and qq are combined and determined from the sideband M bc e (5.21, 5.26) GeV/c 2 , 
AE G (—0.12, 0.12) GeV; we check with MC simulation that the shapes of these backgrounds 
and their ratio in the sideband region are close to those in the signal region. The PDF for 
b— >u is taken from MC simulation. The fraction of p + p~ + pirn is taken from the M hc -AE 
fit; the component f p7T7T alone is taken from the M n ± n0 fit. The fraction oib^u background 
is small and taken from MC simulation. Since fr q q +b _ iC \ = 1 — f pp — f pn7T — fb~*u^ there is only 
one free parameter in the fit. The result is f L = 0.941 ^(mo, wnere the error is statistical. 
Figure 12(b) shows the data and projections of the fit. 

There are eight main sources of systematic error in f L : the pun fraction 

(+0.013,-0.012); the pirn component alone (+0.021,-0.020); the pion identification ef- 
ficiency, which affects the acceptance (+0.000, —0.004); misreconstructed B° — >p + p~ decays 
(+0.005,-0.000); the qq suppression requirement (±0.013); interference of longitudinally 
polarized p's with an S'-wave system in B° — > pun (+0.003, —0.005); a very small bias 

in the fitting procedure measured from a large toy MC sample (+0.000, —0.005); and un- 
certainty in the qq + {b— >c) background shape (+0.004,-0.014). This last uncertainty is 
evaluated by taking the background shape from alternative M bc , AE sidebands. Combining 
all errors in quadrature gives a total systematic error of ±0.030. Thus, 

f l = 0-941 ; ;^ ( s tat) ± 0.030 (syst) . (2) 

To determine CP coefficients A and S, we divide the data into q = ±l tagged subsamples 
and do an unbinned ML fit to their At distributions. Since B 0, s and 5°'s are approximately 
at rest in the T(4S') frame, and the T(45) travels with f3j = 0.425 nearly along the beam 
axis (z), At is determined from the z displacement between the p + p~ and tag-side decay 
vertices: At ps (zqp — Ztag)//?7C- 
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The likelihood function for event % is a sum of terms: 

A = f£! 1>m„ + /& P( At) SCF + /«, P( At) p _ 

<J(J ) 



where the weights are functions of M bc and AE and are normalized to the event frac- 
tions obtained from the M hc -AE and M n±7T o fits. The PDFs V(At) are obtained from MC 
simulation for b— >c and and from an M bc sideband for qq. We include a term for SCF 
events in which a ^ daughter is swapped with a tag-side track; the PDF and function / SCF 
are also obtained from MC simulation. 

The signalPDF is e~ |A * l/r Bo/(4r BO )x{l =p Alo £ ±_(1 - 2u^) [.Acos(Am At) + Ssin(Am At) 
where the upper (lower) sign corresponds to B°(B°) tags, cu e is the mistag fraction for 
the £th bin of tagging parameter r, and Au> e is a possible difference in uj t between B° 
and B° tags. Values of u e and Au;^ are determined from a large P° — > D*~£ + v sample. 
Coefficients A and S receive contributions from longitudinally (L) and transversely (T) 
polarized amplitudes, e.g., A = f L A L + (1 — f L )A T . The transversely polarized ampli- 
tude has a CP-odd component. For a negligible penguin contribution, A T = A L but 

S T=[( l -fh- 2 f cp-odd) I (1 - /l)] ^l; since /cp-odd < /t and /t is small ; we assume A—A L , 
S = S L , and take the possible difference as a systematic error. 

The signal PDF is convolved with the same At resolution function as that used for 
Belle's sin 2<p l measurement The PDFs V p7rn and V SCF are exponential with t = t b and 
r^O.93 ps (from MC simulation), respectively; these are smeared by a common resolution 
function. We determine A and S by maximizing J2i log A; where % runs over the 656 events 
in the M hc -AE signal region that satisfy M w±n0 € (0.62, 0.92) GeV/c 2 . The results are 
A = 0.00 ± 0.30 and S = 0.08 ± 0.41, where the errors are statistical. The correlation 
coefficient is —0.057. These values are consistent with no CP violation (A = S =0); the 
errors are consistent with expectations based on MC simulation. Figure El shows the data 
and projections of the fit. 

The sources of systematic error are listed in Table HJ The error due to wrong-tag fractions 
is evaluated by varying uj t and Aujg values. The effect of a possible asymmetry in b— >c and 
qq is evaluated by adding such an asymmetry to the b — > c and qq At distributions. The 
error due to transverse polarization is obtained by first setting f L equal to its central value 
and varying A T , S T from —1 to +1; then assuming A T = A L , S T = —S L (f T is CP-odd), 
and varying f L by its error. The sum in quadrature of all systematic errors is ±0.09. Thus, 

A L = 0.00 ± 0.30 (stat) ± 0.09 (syst) (3) 
S L = 0.08 ± 0.41 (stat) ± 0.09 (syst) . (4) 



These values are similar to those obtained by BaBar 

We use these values and the branching fractions for B° — > p + p~ ll2Tl. p + p° [HI], and p°p° 0] 
to constrain 2 . We assume isospin symmetry [T^j and follow Ref. |l5l]. neglecting a possible 
1=1 contribution to B° —> p + p~ • ^ e ^ TS ^ ^ ^ e measured values to obtain a minimum 
X 2 (denoted Xmin)j we then scan 2 from 0° to 180°, calculating the difference A% 2 = 
X 2 (02) — Xmin- We insert A% 2 into the cumulative distribution function for the x 2 distribution 
for one degree of freedom to obtain a confidence level (C.L.) for each 2 value. The resulting 
function 1 — C.L. [Fig. E[d)] has more than one peak due to ambiguities that arise when 
solving for (f> 2 - However, only one solution is consistent with the Standard Model 
(88 ± 17)° or 59° < 2 < 115° at 90% C.L. 
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FIG. 3: The At distribution of events in the M bc -AE signal region that satisfy M ± 6 
(—0.62, 0.92) GeV/c 2 , and projections of the fit. (a) q = +1 tags; (b) q = — 1 tags; (c) raw 
CP asymmetry for good tags (0.5<r<1.0); (d) 1 — C.L. vs. (f> 2 - In (a) and (b), the hatched region 
(dashed line) shows signal (background) events. In (d), the solid (dashed) horizontal line denotes 
C.L. = 68.3% (90%). 



In summary, using 253 fb _1 of data we have measured the branching fraction, polarization 
fraction, and CP coefficients A and S for B°-^p + p~ decays, and constrained the angle <f) 2 - 
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program of KOSEF (Korea); KBN (contract No. 2P03B 01324, Poland); MIST (Russia); 
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TABLE I: Systematic errors for CP coefficients A and S. 



Type 


AA ( 


XlO- 2 ) 


AS (XlO- 2 ) 




+ <T 


—rr 


+o- 


—rr 


Wrong tag fractions 


0.5 


0.6 


0.8 


0.8 


Parameters Am, t b0 


0.1 


0.1 


0.9 


0.9 


Resolution function 


1.3 


1.3 


1.3 


1.3 


Background At distributions 


1.6 


1.5 


2.3 


2.5 


Component fractions 


2.1 


2.6 


5.1 


4.5 


Background asymmetry 


0.0 


2.0 


0.0 


4.3 


Possible fitting bias 


0.0 


1.0 


0.7 


0.0 


Vertexing 


4.1 


2.8 


1.3 


1.4 


Tag-side interference [101 


3.7 


3.7 


0.1 


0.1 


Transverse polarization 


6.3 


6.3 


7.1 


5.8 


Total 


+8.9 


-8.8 


+9.3 


-9.2 
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